function [xi1, xi2] = Xi(k, mu, eta, c1, c2, c3)
    % solve_xi 求解超越方程组，计算 xi1 和 xi2
    %
    % 输入参数：
    %   k, mu, eta : 参数
    %   c1, c2, c3 : 常数参数
    %
    % 输出参数：
    %   xi1, xi2 : 方程组的数值解
    
    % 定义超越方程组
    F = @(x) [ 2*x(1) - k*sin(x(1) + mu) + k*sin(mu) - c2 - c3;
               x(1) + x(2) - k*sin(x(2) + eta) + k*sin(eta) + c1 - c2 ];
           
    % 初始猜测
    x0 = [0; 0];
    
    % 调用fsolve求解
    options = optimoptions('fsolve','Display','off');
    sol = fsolve(F, x0, options);
    
    xi1 = sol(1);
    xi2 = sol(2);
end
